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SIMULATIONS OF PROMINENCE FORMATION IN THE 
MAGNETIZED SOLAR CORONA BY CHROMOSPHERIC 

HEATING 
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ABSTRACT 

Starting from a realistically sheared magnetic arcade connecting chromo- 
spheric, transition region to coronal plasma, we simulate the in-situ formation 
and sustained growth of a quiescent prominence in the solar corona. Contrary 
to previous works, our model captures all phases of the prominence formation, 
including the loss of thermal equilibrium, its successive growth in height and 
width to macroscopic dimensions, and the gradual bending of the arched loops 
into dipped loops, as a result of the mass accumulation. Our 2. 5- dimensional, 
fully thermodynamically and magnetohydrodynamically consistent model mimics 
the magnetic topology of normal-polarity prominences above a photospheric neu- 
tral line, and results in a curtain-like prominence above the neutral line through 
which the ultimately dipped magnetic field lines protrude at a finite angle. The 
formation results from concentrated heating in the chromosphere, followed by 
plasma evaporation and later rapid condensation in the corona due to thermal 
instability, as verified by linear instability criteria. Concentrated heating in the 
lower atmosphere evaporates plasma from below to accumulate at the top of 
coronal loops and supply mass to the later prominence constantly. This is the 
first evaporation-condensation model study where we can demonstrate how the 
formed prominence stays in a force balanced state, which can be compared to 
the Kippenhahn-Schliiter type magnetohydrostatic model, all in a finite low-beta 
corona. 
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INTRODUCTION 



Prominences represent fascinating large-scale, cool (~ 8000 K) and dense (10 



10 



10 



ii 



cm -3 ) structures, suspended in the hot and tenuous solar corona above magnetic neu- 
tral lines, which separate opposite polarity photospheric ma gnetic regions. T he magnetic 
field strength in quiescent prominences lies between 3-15 G ( ILeroy et al.lll983l ). It is much 
stronger in active region prominences. The field topology of prominences is less accurately 
known, but was found t o be mostly horizontal, m aking an acute angle with respect to the 
axis of the prominence (jBommier fc Leroyl I1998T 1. It is accepted that the stable existence 
of prominences is due to the mechanical support and thermal shielding of the magnetic 
fields. Two topological l y diff erent types of static promine nce models were proposed by 
Kippenhahn fc Schluterl (119571 ) and lKuperus fc Raadul (I1974I ). In the Kippenhahn-Schliiter 
(K-S) prominences, also known as normal-polarity prominences, the horizontal magnetic 
field through the prominence points from the underlying positive to the negative magnetic 
polarities in the photosphere. In Kuperus-Raadu (K-R) prominences, or inverse-polarity 
prominences, the magnetic field points from negative to positive. In both types, the concave- 
upward parts of magnetic field lines or 'dips' host and support prominence material vi a the 
magnetic tension force against gravity. These have been extended by many authors (ILow 
19751 ; lAmari fc Aly!ll989l ; iHood fc Anzerlll990i ). and recently, flux- rope embedded, normal- 
and i nverse- polarity type equilibrium configuration s have been amenable to numerical mod- 
eling ( Petrie et al.l 120071 ; iBlokland fc Keppenjl201ll ). These models merely consider magne- 
tohydrostatic force balance arguments, and leave out the thermodynamic aspects. 

The most elusive aspect in prominence physics is their seemingly in-situ formation. 
One of the earliest suggested formation scenarios relates the sudden appeara nce of promi- 



nence mate rial to a radiative condensation process due to thermal instability ( jParkerlll953 



Field! Il965l ). In the optically thin solar corona, the radiative losses are proportional to 
density squared and the temperature can evolve into the cri tical range where a decrease 
of temperature leads to an increase of the radiative losses (IColgan et al.l 120081 ) . When 
the plasma is perturbed to become dense or cool enough for radiative losses to domi- 
nate both the heating by thermal conduction and any background heating process, sudden 
strong cooling can cause sufficient pressure gradient to suck in material around along field 
lines. Hence, plasma condenses into cool regions and settles in magnetic dips to form large- 
scale prominences. The linear 'trigger' of such thermal- instability condensations was stud- 



ied in multi-dimensional sh eared magnetic field by many authors (Ivan Hoven &; Moklll984 



van der Linden et al.lll992l ). The nonlin ear evolu t ion o f such thermal instability can only 
be researched by numerical simulations. iHildnerl (119741 ) firstly performed two-dimensional 
(2D) numerical simulations of the thermal instability for the formation of pro minences in a 



uniform magnetic and gravitational field, while neglecting thermal conduction. ISparks et al. 
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( 1l990h simulated a nonlinear condensation in a force-fre e sheared magnetic field including 
anisotropic thermal conduction, but neglected gravity. iDrake et al.l ( 119931 ) simulated the 
formation and stable support of a prominence at the apex of a potential magnetic arcade 
including the upper chromosphere. But in their simulations, the initial trigger perturbations 
were artificially added without any detailed mechanism ultimately giving rise to condensation 
into a prominence. 

Many mechanisms for the transport of plasma from the chromosphere to the corona 
thereby triggering the form a tion of prominen ces have meanwhile b een proposed(see the re- 
view by lMackay et al.ll2010f ). lAn et al.l ( 119881 ) and lWu et al.l (119901 ) simulated the formation 
of a prominence in a potential field arcade in an initially i sothermal co r ona b y injecting 
high density material ballistically from the bottom boundary. IChoe fc Led (119921 ) performed 
a 2.5-dimensional (2.5D) magnetohydrodynamic (MHD) simulation to investigate promi- 
nence formation due to siphon flows by photospheric shearing motions. These simulations 
did not include photosphere or chromosphere regions. One-dimensional (ID) evaporation- 
condensation model, which include c hromosphere and chromospheric heating, have been 



studi e d numerically b y many authors (jMok et al.lll990t lAntiochos et al.lll999l ; Karpen et al. 



200ll ; Kia et al.ll201ll ). This thermo-hydrodynamic model depends on heating localized in 
the chromospheric feet of a pre-shaped loop, which evaporates chromospheric material into 
the corona. The increase in density results in a dominating radiative cooling, and the ther- 
mal instability trigger s condensations in the corona, eventually forming a steadily growing or 
recycling prominence (IXia et al.l 1201 if ). In this Letter, we present the first simulation where 
this model is performed in a 2.5D full MHD regime in a realistically stratified sheared arcade 
going from the low chromosphere to high corona. We show how we can simulate the rapid 
thermodynamical changes in a scenario where we form and reach a stable support for a K-S 
quiescent prominence. 



2. NUMERICAL SETUP 

We perform our simulation in a rectangular, Cartesian geometry, ignoring the curva- 
ture of the solar surface. Since prominences are often observed in a horizontally elongated 
form along the magnetic neutral line, we ignore the variation of physical variables along 
the prominence axis (^-direction) and consider their dependence on horizontal x-coordinate 
(perpendicular to the magnetic neutral line) and vertical ^-coordinate, while keeping the 
z-component of any vector. Therefore, we use a 2.5D thermodynamic MHD model including 
gravity, field-aligned heat conduction, radiative cooling and heating terms. The governing 
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equations are as follows 

^ + V-(pv) = 0, (1) 
^1 + V ■ (pvv + p tot I - BB) = pg, (2) 

dE 

— + V • (£v + Pw v — BB • v) = pg • v + V • (k • VT) — Q + H, (3) 

_ + V .(vB-Bv) = 0, (4) 

where p, T, v, B, and I are the plasma density, temperature, velocity, magnetic field, and 
unit tensor, respectively; p to t = p + B 2 /2 is the total pressure, composed of thermal pressure 
p and magnetic pressure B 2 /2; E = pj (7 — l)+pv 2 /2+B 2 /2 is the total energy density, where 
7 = 5/3 is the ratio of specific heats; g = — g -R©/ (R Q + y) 2 y is the gravitational acceleration 
with Rq the solar radius and go the solar surface gravitational acceleration; n is the thermal 
conductivity tensor; Q and H are the radiative loss rate and heating densities, respectively 
Considering fully ionized plasma with 10:1 abundance of hydrogen and helium, we have 
p = 1.4m p nH, where m p is the proton mass and nn is the number density of hydrogen. 
As for the equation of state, we adopt the ideal gas law p = 2.3n il k B T, where k B is the 
Boltzmann constant. The radiative cooling term is taken as Q = 1.2n^A(T), where A(T) is 
the radi ative loss functi on for optically thin emission, which was also used in our earlier ID 



studies (IXia et al.ll201ll ). The anisotropic thermal conduction along the magnetic field lines 
is included with k = /tbb, where b is the unit vector along B and k = 10~ 6 T 5//2 erg cm -1 
s _1 K _1 is the Spitzer conductivity. 

As for the initial magnetic field, we intend to mimic sheared arcades above a neutral 
line {x = 0, y = 0) and start our simulation from a force-ba lance d state . Therefore we adopt 



an analytic solution of a nonlinear force- free field found by iLowl (119771 ) as follows 



B x = 2B ( ky + Lj^j //, B y = -2B kx/f, B z = ^^/f, (5) 



where f= M £ + k 2 x 2 +(ky + — ^ ) , (6) 



,7^ + k 2 x 2 + (ky + \— ^ 

where |p| < 1 controls the shearing rate of the arcades (no shearing if p = 0) and k controls 
the spatial concentration of the field. The shearing decays from lower loops to higher loops. 
We set p = 0.95, k = 0.5, and B = 4 G, leading to a realistic 2.5D arcade topology, where 
the field lines at the height of 20 Mm make an angle of 45° with the underlying neutral line. 



For the initial thermal structure, we set a chromosphere with temperature of 10000 
K below a height of 2.7 Mm and choose a vertically stratified temperature profile with a 
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constant ve rtical thermal conduction flux (i.e. n dT/dy = 2 x 10 5 ergs cm 2 s x ) above 
the height (IMok et al.l l2005t iFontenla et al.lll99ll ). The initial density is then determined 



by assuming a hydrostatic atmosphere with the number density of 2 x 10 cm at the 
bottom. To use a proper background heating term to maintain a hot corona, we are inspired 



by pa rametric comparisons of different models of coronal he ating done by iMandrini et al. 



( 120001 ) and assume the heating rate to be proportional to B 2 (IMok et al.ll2008l ). However, in 
order to compensate the radiative loss in the transition region above the neutral line, where 
dominating horizontal field lines insulate this region from getting heat thermally conducted 
from the corona above, we add an extra heating equal to the local radiative loss, below a 
6 Mm height, purely concentrated in this region. The resulting two- component, parametric 
background heating is expressed as 



H = c B 2 + 0.5Qcos 



7rR 



tanh 



with c = 7.5 x 10 



6 erg cm 3 s 1 G 2 , 



Ilk 



2B 

= 6 Mm, and yd 



y-yh 

Vd 

0.5 Mm. 



(7) 



The configuration of this system is symmetric about the y-axis. We exploit the symme- 
try to study prominence formation under symmetric heating conditions, allowing to reduce 
computational domain to the right half of the simulated area (within < x < 30 Mm and 
< y < 40 Mm). Symmetric/asymmetric boundary treatments can then be used at the 
y-axis, while for right-hand side boundaries we adopt a zero velocity, continuous density and 
pressure, and fixed magnetic field. The top and the bottom boundaries have a zero velocity, 
fixed magnetic field and extrapolated density and pressure, respectively, which are derived 
by assuming a hydrostatic equilibrium with the temperature being fixed. 

We use the parallelized Ad aptive Mesh Refinement Versatile Advection Code (MPI- 
AMRVAC) (IKeppens et al.ll2012l ) to solve the governing equations with a second-order shock- 
capturing Total Variation Diminishing Lax-Friedrichs scheme. The effective resolution of 
512 x 1024 is attained by using 5 levels of AMR. The equivalent spatial resolution is then 
39 km/59 km in the vertical/horizontal direction. Anisotropic thermal conduction is added 
explicitly as an energy source term. 

This initial state is not in thermal equilibrium, and we integrate the governing equations 
in time with the background heating H = Hq until the system relaxes to a quasi-equilibrium 
shown by panels (a) and (b) in Figure [H The chromosphere of about 4 Mm thickness is 
connected to the corona by a very thin transition region. The projected magnetic field lines, 
colored according to the local density, are plotted through selected footpoints at the bottom 
and side boundaries. The plasma beta is 0.1 at 20 Mm height above the neutral line while 
the temperature and number density are 1.6 MK and 2.6 x 10 8 cm -3 there, respectively. The 
maximal residual velocity is small, less than 5 km s" 1 . Starting from this quasi-equilibrium, a 
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relatively strong heating Hi is added. This extra heating is localized near the chromosphere 
(see the contours in Figure (Ha)) with its formula as: 

H = iC^By/Bof ify^yc, 
' [C 1 ( J B y / J B ) 2 exp(-(y- 2/c ) 2 /A) ify>y e , 

where C% = 10 _2 ergs cm -3 s -1 , y c = 3 Mm, and A = 3 Mm. This localized heating is 
concentrated in the regions of strong B y . 



3. PROMINENCE FORMATION DUE TO EVAPORATION AND 

THERMAL INSTABILITY 

As the localized heating is functioning, chromospheric plasma is heated and evaporated 
into the arched coronal loops, increasing the density and the temperature there. About 6 
minutes later, the temperature reaches a maximum value of 2.2 MK, and then starts to 
decrease slowly (see Figure [T](c, d)). At about 84 minutes, the temperature of an inverted 
triangle-shaped region around the apexes of a bundle of magnetic loops decreases drastically 
(see Figure[T](f)) and a small condensation with the typical chromospheric density (6.3 x 10 10 
cm -3 ) appears near the loop top at a height of 25.4 Mm. This is accompanied by two strong 
inflows moving towards the central condensation segment, with a maximum velocity of 70 km 
s _1 from the two sides (see Figure [He)). The magnetic field lines near the condensation make 
an angle of 62° with respect to the ^-direction. After this sudden birth of the prominence, 
condensation successively happens on the tops of lower and higher coronal loops, leading to 
a rapid extension of the prominence in the vertical direction. As the inflows collide near the 
loop tops, two rebound shock waves are formed (see diamond-shaped wave fronts in Figure 
□Jg, h)) and propagate from the apex towards the loop feet, during which they sweep across 
and slow down the evaporated upflows. 

In order to investigate the thermal instability during the in-situ formation of the promi- 
nence, we quantify the temperature, the time derivati ye of temper ature, the density, the 



pres sure, the th ermal instability isochoric criterion Cp (jParkerl Il953l ) and isobaric criterion 



C F fjFieldl lliiih at the site (x—0, y=25A Mm) of the first condensation (see Figure [2]). A 



value of 20 Mm is adopted as the wavelength of perturbation when calculating the thermal 



instability criteria (IXia et al.ll201ll ). After £=82.3 minutes, when the temperature and pres- 
sure start to decrease nonlinearly, both Cp and Cp dive into the negative region indicating 
the functioning of thermal instability. Although the localized heating in the chromosphere 
is in a wide range in the x-direction, plasma condensation appears only on long field loops. 
The plasma conditions along these field loops satisfy the criterion of thermal instability. At a 
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given height, a larger shear rate (quantified by /i) leads to a longer field line whose footpoints 
are closer to the neutral line, which renders the sheared arcade more vulnerable to thermal 
instability. 



4. OVERALL PROMINENCE STRUCTURE AND FORCE BALANCE 

As the prominence grows fatter and heavier, the field lines penetrating the prominence 
are gradually bend downwards forming vertically aligned dips. This is shown in Figure [3f a), 
a 3-dimensional (3D) illustration of the vertical 'sheet-like' prominence and selected field 
lines at t — 143 minutes. The right panel of Figure [3] gives the corresponding 2D projected 
image. At this moment, the prominence has a vertical extension of 12 Mm from its bottom 
at 15.3 Mm to its top at 27.3 Mm. The horizontal width of the prominence is about 2.7 Mm 
at the middle and very thin at top and bottom edges. In the prominence, the temperature 
is about 18600 K, the density varies from 7.3 x 10 10 at the bottom to 4.3 x 10 10 cm at the 
top, and the plasma beta changes from 0.22 to 0.47 as the magnetic field strength changes 
from 6.8 G at the bottom to 3.7 G at the top. The coronal loops below the prominence 
are heated to nearly 2 MK, which is hotter than the other coronal area (~1.4 MK). The 
transition region a bove the neutral line has risen to a hei ght of about 7 Mm, which is similar 



to previous works (jLionello et al.ll2001t iMok et al.ll2005l ). 



The total mass of the prominence is increasing as the localized chromospheric heating 
is kept. In order to quantify the prominence mass, we integrate the plasma with density 
larger than 4.2 x 10 10 cm -3 within a square box which surrounds the prominence. Within 
the first 10 minutes, the prominence mass grows nonlinearly with a mean rate of about 764 
g cm -1 min -1 . Later the prominence mass grows linearly, with a rate of about 519 g cm -1 
min -1 . The temporal evolution of the total mass of the prominence in a unit length in the 
^-direction is quantified in Figure HI 

The initial magnetic field was force-free, with the current being parallel to the mag- 
netic field. After the formation of the prominence, the current is locally increased in the 
prominence region. The ^-component of the current density in the prominence is positive 
and is significantly stronger than in the surrounding corona, as shown by Figure [3(a). The 
magnetic field in the plane is pointing from the left to the right, so the Lorentz force in the 
prominence is pointing upward and is able to balance the gravity of local dense plasma. This 
is convincingly demonstrated by comparing the gravity, Lorentz force, gas pressure gradi- 
ent, and their sum, through the prominence structure. We check this first along a central 
vertical line si through the prominence (see Figure Ufa)). The distributions of these forces 
are displayed in Figure Efb). Near the prominence-corona transition regions (PCTR) at 15.3 
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Mm and 27.6 Mm heights, the gas pressure gradient and Lorentz force fluctuate rapidly, as 
this transition region is complicated by the local thermodynamics driven by the radiative 
losses combined with field-aligned heat conduction. In this thin PCTR layer, forces are 
not in balance exactly, but the average value of the resultant force to remain close to zero. 
However, in the prominence body, the dominant Lorentz force and gravity nearly balance 
perfectly, except for small fluctuations of the Lorentz force, which are then compensated 
by opposite fluctuations of the pressure gradient. The resultant force vanishes throughout 
the prominence body, as indicated by the solid line, realizing force-balance in the vertical 
direction. In the horizontal direction, the distributions of the Lorentz force, the pressure 
gradient, and their resultant force along the slice s2 are plotted in panel (c) of Figure |5j 
We find that inside the prominence (0 < x < 1.1 Mm), the force of the pressure gradient 
points outwards and almost balanced by the Lorentz force, which is pointing to the center 
of the prominence. Hence the horizontal balance is realized by a magnetic pinching. The 
strength of these forces increases from the center to the edge of the prominence, and a small 
resultant force points outwards. In the thin PCTR region, forces fluctuate again indicating 
the detailed thermodynamic processes at play there. In the end, the mass of the prominence 
would be saturated, after which any newly-formed condensation would fall aside from the 
corresponding magnetic dips, forming the drainage of cool material. 



5. CONCLUSION 

In this Letter, we simulate the formation of a normal-polarity quiescent prominence 
in a magnetized coronal arcade by chromospheric heating for the first time in a realistic 
multi-dimensional magnetic configuration. There should be no fundamental difference in the 
inverse-polarity configuration, which requires 3D simulations, as far as the radiatively driven 
condensation is concerned. The magnetic dips supporting the prominence mass against grav- 
ity are self-consistently formed in an overall low-beta environment. Therefore, the magnetic 
dips in our model are a consequence of the prominence formation. Our simulation captures 
many details relevant to quiescent prominence models, and closely resembles many observa- 
tional features. The prominence body, elongated along the invariant direction, is situated 
above the magnetic neutral line and the magnetic field through the prominence makes a finite 
angle to its axis. Our model naturally produces a vertical 'sheet-like' prominence rather than 
multiple threads aligned along the magnetic field lines. As long as the chromospheric heating 
is active, the prominence grows in horizontal and vertical size, while the mass-loaded arcade 
field loops realize a force balance between Lorentz force and gravity throughout the promi- 
nence body. Our simulation radically improves all earlier evaporation-condensation studies, 
which assume rigid one- dimensional field line shapes along which catastrophic cooling see- 
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narios have been studied in parametric detail, invariably leading to field-aligned prominence 
threads. Similar parametric studies will be required in our fully multi-dimensional settings, 
to determine their impact on the macroscopic parameters like prominence width, height and 
total mass. 

This study can guide future observations dedicated to find direct evidence of plasma 
condensation and mass-supply mechanism in prominences and ultimately uncover these mys- 
teries of prominences. Our model can act as a starting point for future studies of Rayleigh- 
Taylor instability d evelopment in quiescent prominences, extending recent local box studies 



( IHillier et al.ll201l[ ). by allowing for true sheared field configurations. We intend to con- 
trast synthetic views on the obtained prominence structures with modern observations, and 
will investigate how spatio-temporally varying heating conditions may give rise to multiple 
condensation sites and filament threads. 
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Fig. 1. — Snapshots of the formation process at t — (first row), 54.4 mins (second row), 
84.4 mms(third row), and 87.3 mins (bottom row). In the left column, the density is shown 
in gray and the projected velocity field is shown by arrows. The right column shows the 
temperature and the projected field lines are colored by the local density. Contours of the 
localized heating are plotted in panel (a) in dotted lines. 



- 13 - 




80.5 81.0 81.5 82.0 82.5 83.0 83.5 84.0 



t (min) 

Fig. 2. — Temporal evolutions of the temperature, the density, the pressure, the thermal 
instability isochoric criterion Cp, and the isobaric criterion Cp at the first condensation site 
(solid lines), as well as the time derivative of the temperature (dashed line, top panel). The 
horizontal dashed lines in the last two panels denote the zero value. 
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Fig. 3. — A K-S prominence at t — 143 mins. (a) 3D impression of the prominence shown 
by the number density and selected magnetic field lines. The bar above the color bar shows 
the opacity of corresponding colors, (b) the temperature and the projected field lines colored 
by the local density. 
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Fig. 4. — Temporal evolution of the prominence mass. 
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Fig. 5. — Panel (a): Distribution of ^-component of the current density (j z ); Panel (b): 
Distributions of various forces along the vertical slice si in panel (a), including the gas 
pressure gradient (dotted line), Lorentz force (dashed line), gravity (dashed dotted line), and 
their sum (solid line); Panel (c): Distributions of various forces along the horizontal slice s2 
in panel (a), including pressure gradient (dotted line), Lorentz force (dashed line), and their 
sum (solid line). 



